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Abstract 

The anomalous properties of the Hall constant in the normal state of high- 
Tc superconductors are investigated within the single-band Hubbard model. 
We argue that the Mori theory is the appropriate formalism to address the 
Hall constant, since it aims directly at resistivities rather than conductivities. 
More specifically, the frequency dependent Hall constant decomposes into its 
infinite frequency limit and a memory function contribution. As a first step, 
both terms are calculated perturbatively in U and on an infinite dimensional 
lattice, where U is the correlation strength. If we allow U to be of the order of 
twice the bare band width, the memory function contribution causes the Hall 
constant to change sign as a function of doping and to decrease as a function 
of temperature. 



PACS numbers: 
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I. INTRODUCTION 



Since the discovery of high-Tc superconductors ten years ago, the anomalous properties 
of their normal state have been the subject of intensive theoretical work. It is widely 
believed that a model of strongly correlated electrons already captures the basic ingredients 
of the relevant physics. In these models, the correlations are represented by a strong local 
interaction U. However, a coherent description of all anomalous properties on the basis of 
such a model is still lacking. The main problem is that exact calculations are generally 
feasible only in a small parameter regime and that most approximation schemes fail in 
capturing all aspects which are supposed to be important. 

The Hall constant is especially hard to describe. One reason for this is that the Hall 
conductivity contains a three-point correlation function after it has been expanded to first 
order in the magnetic field. Then, the calculation of vertex corrections is a tough problem 
which, to our knowledge, has been attempted only in the case of a Fermi liquid and to 
leading order in the quasiparticle damping |l|]. Moreover, since the frequency dependent 
Hall constant is given as a quotient of conductivities, the limit — may be precarious due 
to resonances like the Drude peak. A more technical peculiarity of the Hall effect is due to 
the fact that the magnetic field is introduced via a vector potential which, formally, breaks 
the symmetry with respect to lattice translations. But even in the simplest case of a Bloch- 
Boltzmann description, the temperature dependence of the Hall constant may be difficult 
to reproduce, because the relaxation time cancels once it is assumed to be independent of 
momentum. 

The measurements of the Hall constant in high-T^ materials reveal two major anoma- 
lous dependences: on temperature and on doping. Both cannot be understood within con- 
ventional band theory. For noninteracting tight-binding electrons on a two-dimensional 
square lattice, the Hall constant changes sign at half filling as the Fermi surface changes 
its shape from electronlike to holelike. In contrast, the Hubbard model in the large-?/ limit 
exhibits an additional sign change below half filling which is purely due to correlations . 
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In addition, in the limit 5 — > 0, i.e. near half filling, the Hall constant diverges according 
to a 1/5 law 0,^. These properties are supposed to account for the doping dependence 
observed in, e.g., La2-xSrj:Cu04 As for the anomalous temperature dependence of the 
Hall constant, the most striking features are: firstly, a strong decrease which is, in some 
cases, as fast as 1/T 0; secondly, the lack of saturation above a fraction of the Debye 
temperature, typically ~ 0.2 — 0.4Tci in contrast to what is expected in a Fermi liquid 
description with weak electron-phonon coupling 0; and thirdly, a quadratic dependence of 
the Hall angle on temperature for not too large dopings [l8|, [T0| , [TT| . In a Fermi liquid, the 
temperature dependence arises from an anisotropic relaxation time . If we assume scat- 
tering off phonons to be the main inelastic process, a temperature dependence is conceivable 
only below a certain temperature scale: Then, a sufficiently anisotropic Fermi surface causes 
the scattering to be confined to those regions of the Fermi surface, where small momentum 
transfers are possible. For high enough temperatures, this kinematic restriction is lifted and 
the scattering becomes isotropic, thus leading to a cancellation of the relaxation time. The 
cross-over temperature is given by ~ 0.2 — 0.4Td. The universally observed decrease of the 
Hall constant as a function of temperature in almost all high-T^ compounds up to tempera- 
tures clearly beyond this temperature scale must therefore be due to electronic correlations 
as well. 

In the following, we investigate the Hall effect on the basis of the simplest model of 
strongly correlated electrons, namely the single-band Hubbard model on a hypercubic lattice 
in d dimensions with nearest-neighbour hopping. This model, along with Mori's formalism 
used to represent the Hall constant, is introduced in Sec. II. In this theory, the Hall constant 
is given as the sum of its infinite frequency limit {RJj) and a memory function contribution. 
The former term was first considered by Shastry, Shraiman and Singh 0. Our emphasis is on 
the memory function term which represents the deviation of the Hall constant from for 
finite frequencies and thus cannot be neglected when considering the case of zero frequency. 
The advantage of our representation of the Hall constant is that we do not have to cope 
with a quotient of conductivities as opposed to the usual approaches. This is why the Hall 
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constant at low frequencies becomes less sensitive to the detailed resonance structure of the 
conductivities. In the remainder of this paper, this advantage is exploited for the range of 
weak to intermediate correlation strengths, while the opposite limit of strong correlations will 
be addressed in a forthcoming paper ||T2|. In Sec. Ill, we proceed by calculating the memory 
function to second order in the Hubbard interaction and to first order in the magnetic 
field. Expansion with respect to the magnetic field leads to a decomposition of the memory 
function into two terms, namely a two-point and a three-point correlation function. Both 
contributions are evaluated exactly in infinite spatial dimensions. Our results indicate that 
the memory function term is important. Only then, a precursor effect of the sign change of 
the Hall constant as a function of doping appears even in perturbation theory. Moreover, 
when extrapolating our results to U values of the order of twice the free band width W, we 
get most of the qualitative features observed in, e.g., La2_xSr2;Cu04: the sign change with 
respect to doping and the decrease of the Hall constant up to unusually high temperatures, 
characteristic of most high-T^ compounds. Finally, in Sec. IV, we summarize our main 
results. 



The single-band Hubbard model on a d-dimensional hypercubic lattice in a magnetic 
field reads 



II. THEORETICAL FRAMEWORK 



A. Single-band Hubbard model 



H 



T + V, 



(1) 



T 



(2) 



<i3> 



V 



(3) 



where the sum in the hopping term T is restricted to nearest neighbours and V is the 



Hubbard repulsion. The Peierls phase factor 
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gauge invariance [0 and the sign of the charge e is chosen to be negative. Since only 



nearest neighbour hops are taken into account, we may approximate 

where Ri denotes the lattice vector to site i. The vector potential decomposes into two terms 
describing the electric and magnetic field, respectively: 

A{t, f) = A^^ (t) + A"""^ (f ) , (5) 
m = -^/''\t) , (6) 
H = rotA""''{f) . (7) 

In linear response theory with respect to the electric field, the latter appears only in the 
definition of the current operator. More precisely, the current operator is defined as the 
following functional derivative: 



^ _1 6H{t) 



(8) 



The homogeneous magnetic field is chosen to point in the z-direction and it is advantageous 
to fix the gauge from the very beginning according to the Landau choice 

A-^-3[R) = R^Ey^ (9) 

since then, translational symmetry is broken only in one dimension, namely the x-direction. 
y is a primitive lattice vector. We need the current operator only up to first order in the 
magnetic field: 

J, = Ji'^^ + sX , (10) 

J2 S,cUR + 6)c,{R), (11) 

R,R+5,cT 

6l = -et J2 SjA^''>^{R)ct{R + S)c„{R) . (12) 

R,R+5,c7 

Here, 5 is a nearest neighbour vector and the summation is over pairs of nearest neighbours. 
Note, however, that due to the gauge fixation (|), we cannot choose periodic boundary 



conditions in the x-direction. Thus, if S points in the x-direction, we have to carry out the 
sums in such a way that the components and R^+S^ are simultaneously elements of the set 
consisting of the x-coordinates of all lattice sites, e.g. {-R^™, -R^™ + 1, . . . , -R™"^}. Here, it is 
assumed that the lattice has N^j. sites in the x-direction which implies R^"'^ = R^^^ + N^ — 1. 
Of course, observable quantities are not allowed to depend on the lattice location i?^*"^. The 
hopping term is expanded analogously, yielding: 

f = f + 5f , (13) 
fW = -t Y: ctiR + 5)c.{R), (14) 

R,R+5,a 

5f = -iet SA"''''{R)ct{R + S)c^{R) ■ (15) 

R,R+5,cr 

The term without magnetic field, i.e. Eq. (|T^), becomes diagonal in crystal momentum space 
with a band dispersion = —2tY,'j=i coskj. 



B. Mori theory 

In this subsection, the basics of Mori's memory function formalism is reviewed briefly. 



For further details, see e.g. Ref. |T^. The best known application of Mori theory is the 



description of many particle systems in the hydrodynamic regime |jT5|. There, one is only 



interested in the dynamics of the hydrodynamic variables. They are characterized by the 
fact that their transport is restricted by conservation laws or by broken symmetries. Thus, 
they are bound to vary on a time scale that is very slow in comparison to that of all the 
other degrees of freedom. Now, the Mori theory enables one to separate these two time 
scales: The equations of motion of the hydrodynamic variables take on the form of coupled 
integro-differential equations. The corresponding integral kernels of these so called Mori 
equations are memory functions in which the influence of all the other degrees of freedom 
is accumulated, hence the name "memory function" . In this context of hydrodynamics, the 
memory functions are rapidly varying functions, whose effect may be simulated by damping 
constants. Then, the Mori equations take on a form analogous to that of the Langevin 

6 



equation for a particle undergoing Brownian motion. However, the validity of the Mori 
equations is not restricted to the special set of hydrodynamic variables. In the simplest case, 
one sets up the Mori theory for those observables that constitute the correlation functions 
one is interested in. This leads to representations for the unknown correlation functions in 
terms of memory functions in which all analj^ic properties are fulfilled by construction. On 
the other hand, it may be difficult to find an approximate expression for a given memory 
function. 



The Liouville space C is defined as the linear vector space over the field of complex 
numbers whose elements are the linear operators in the familiar Hilbert space of quantum 
mechanics, and where the usual operations like scalar multiplication etc. hold. In this 
Liouville space exist linear operators that are called superoperators to distinguish them from 
the usual ones. (Henceforth, normal operators are denoted with a hat, superoperators not.) 
The most important superoperator is the Liouville operator L, which maps a given operator 
onto its commutator with the Hamiltonian: 



An other important class of superoperators are superprojectors. However, their definition 
implies a scalar product in C In the context of response functions, the most convenient 
scalar product turns out to be the so called Mori product 



where < . . . > denotes the thermal average and /3 is the inverse temperature. On the basis 
of this scalar product, we may now speak of adjoint superoperators S and S~^, and thus 
of unitary and Hermitian ones in the usual sense. The projector P that projects onto the 



1. basic notions 



LA = [H, A] . 



(16) 




(17) 



subspace of C spanned by linearly independent elements \Gi), reads: 
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P = Y.\G^)g^AG,l (18) 

where the metric gij is the inverse of the matrix [Gi\Gj)^ i.e. X^fc 9iki,Gk\Gj) = Sij. In fact, 
this imphes the idempotence property = P. Finally, the definition of the Mori product 
implies the validity of the so called Kubo identity 

(3iA\L\B)=< >, (19) 

which will play an important role. 



2. Memory function approach to the Hall constant 

From Eq. (|19D follows a representation for the current- current correlation function <^ 
Ju J^i , defined as the Laplace transform of —i < [d^^Jy , J^] >: 

^X-Jt. >.= -l3{JA^^\Ju) = -X^uiz) . (20) 

Here, z is a complex frequency, which ultimately has to be specialized to u + iO^. For 
formal manipulations, however, it is more convenient to deal with the complex frequency z 
rather than with u. The last expression has to be inserted into the Kubo formula for the 
conductivity tensor, 

"^"^^^^ = ^{< r^f.> + ^Ju; Jf, > J , (21) 

where < f^j^ > arises from the equilibrium part of the current and is defined as the second 
functional derivative of the Hamiltonian with respect to the external electric field: 

1 6^H{t) 



(22) 



e2 6Afj{t) M^'(t) 

In the Hubbard model its expectation value is given as < f^fj, >= ItY,]:^ 

cos kx < 

ng^ >, i.e. as the average kinetic energy per dimension. Using the fact \imz_,Q z (T^^{z) = 
which holds for a metal in the normal state, we show that < f^^ > also equals the static 
susceptibility defined through XfMu{z = 0) = x° • 
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X° = =< >=-^<f> . (23) 

Thus, the conductivity tensor may be written as 

a,,{z) = ^(3{J^\-^\l) ^ ^C,,{z) . (24) 

In order to represent the relaxation functions C^^{z) in terms of memory functions, we 
introduce the superprojector P that projects onto the subspace of £ spanned by the current 
operators Jx and Jy-. 

X u=x,y 

and the complementary superprojector Q = 1 — P. By making use of the operator identity 
^-T = - — -6— rr with a = z + LQ and b = LP, we find: 

a+b a a a+b ^ ' 

Cf^uiz) = - lR,a{z)C^Az) , (26) 



R.Az)^^,i.U^L\l) (27) 



= + iM^,{z) . (28) 
The terms in the last equation are the frequency and the memory matrix, respectively: 

n^, = ^< [J^ ,Ju]>, (29) 

M,Az) ^ ^iQLJ,\ 'j QLl) . (30) 

The last equation shows that the dynamics of the memory functions is governed by QLQ 
rather than L. Solving Eq. (^Gf ) for the matrix C{z) leads to 



C{z) = ix\zl + Vt + i'M{z))-\ (31) 

Together with Eq. (|2^) , this demonstrates that the Mori theory enables us to calculate 
directly the resistivity tensor. Therefore, the desired representation for the dynamical Hall 
constant can be read off from the last equation: 

Rh{z) = -^riT 1™ — — TT ■ (32) 
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Since the memory function term will be shown to vanish in the high frequency limit as 1/^;^, 
the first term represents the high frequency limit of the Hall constant considered by Shastry 
et al. 0]: 

i?- = hm %^ . (33) 

Moreover, is the generalization of the cyclotron frequency to the lattice case. Therefore, 
within a Boltzmann equation approach, only the term (|33D is considered. The goal of the 
subsequent sections is to investigate the memory function term Mi^z) = Mn^yi^z) for finite 
frequencies. 

3. Analytic properties 



The analytic properties of M(z) in Eq. ( p2D may all be derived on the basis of Eq. (pOl). 
An alternative procedure is to solve Eq. (|31|) for NL{z) and go back to the analytic properties 
of the current susceptibilities Xfj,u{z), cf. Eq. pPj). M{z) reads in terms of the susceptibilities 

'm^) = , - ^Mii . (34) 

From time reversal invariance, homogenity of time and the fact that the current operators 
are Hermitian, we may deduce the following symmetry properties ||15||: 



Xxx{-z) = Xxxiz) , (35) 

Xlxiz) = Xxx{z*) , (36) 

Xxy{-Z) = -Xxy{z) , (37) 

Xlyiz) = -Xxyiz*) . (38) 

Together with Eq. (^), this implies: 

M{-z) = M{z) , (39) 

M*{z) = M{z*). (40) 
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M{z) can be represented spectral integral 

rd.m^^ (41) 

J TT UJ — Z 

where the spectral function M"{uj) is given by the discontinuity across the real axis: 

M{uj ± iO+) = M'{uo) ± iM"{uj) . (42) 

From the analytic properties (^OD, it follows that M'{uj) and M"{uj) are real functions sat- 
isfying 

M'{-uj) = M'{uj) , (43) 

M"{-uj) = -M"{uj) . (44) 

Thus, two further conclusions can be drawn: Firstly, only even powers m 1/z contribute 
to the high frequency expansion of M{z). And secondly, the quotient M"{uj)/uj must be 
integrable around u; = 0, 

dujM"{uj)/u <oo, (45) 

-J 

which can be seen from the fact that the dc-Hall constant contains this integral. Note, that 
we need not understand this expression as a Principal value integral due to the fact that the 
integrand is even. 



III. PERTURBATION THEORY 

Despite many interesting works on the normal state Hall effect of high-Tc superconduc- 
tors, a calculation that incorporates all the complicated many-body correlations exactly 
within a microscopic model is still lacking. The following treatment of the Hall constant 
closes this gap at least in the perturbation-theoretical regime. On the other hand, the rel- 
evant parameter regime is believed to be the strong correlation limit rather than the weak 
one. However, it turns out that the final expression may well describe the observed depen- 
dencies on temperature and on doping at least qualitatively, if we allow U to be extrapolated 
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to values of the order of the free band width W. Thus, a precursor effect of the anomalous 
dependences clearly shows up even in the regime of weak correlations. 



A. Approximation 

The perturbation-theoretical treatment of the Hall constant is by no means straightfor- 
ward. As is well known, the evaluation of response functions like Xfiu{z) of Eq. (|20| ) by 
expansion in a small interaction parameter fails because, as an artefact of such an expan- 
sion, these functions become singular for small frequencies z. This difficulty was resolved 



by Gotze and Wolfle |jT6[ some time ago by means of a memory function approach. They 
calculated the memory function perturbatively, which, at first, is valid only at high enough 
frequencies. It turns out, however, that their expression for the memory function depends 
only smoothly on frequency below a certain frequency scale and tends to a constant in the 
limit uj 0. Furthermore, in their approximation scheme the correct resonance structure 
of the studied response functions is inherently built in. Thus, their results could be used 
in the whole frequency regime including the hydrodynamic one. However, we cannot carry 
over their analysis straightforwardly to the present problem, because otherwise, we would 
encounter a spurious singularity in the limit that is of most interest, i.e. — > 0. 

In this subsection, we identify the precondition which is necessary to obtain regular ex- 



pressions in this limit and that was fulfilled trivially in the applications of Ref. [0 but 
is not in our case. Since this condition does not affect the correct description of the local 
interaction U, we may take it as an approximation. By using Mori's formalism, we shall 
see that once this condition is assumed to be satisfied no further approximations have to 
be made. This last point cannot be seen in the more intuitive introduction of the memory 



function concept as given in Ref. [jl6[ and shows that the extrapolation to low frequencies 
therein is exact. 

Perturbation theory is based on the following decomposition of the Liouville operator: 
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L = Lo + L,, (46) 

where Lq and Li are assigned to the hopping and interaction term of the Hubbard hamilto- 
nian ([1|), respectively. The perturbation-theoretical regime is given by the condition U <^ W, 
where W is the width of the free band and thus represents the characteristic energy scale 
introduced by Lq. The precondition to obtain a regular expression for the memory function 
M{z) for all frequencies to leading order in U is that the relevant operators Jx and Jy span 
a subspace of C that is invariant with respect to actions of Lq [|l4l. If this condition were 
satisfied in our case, it would take on the form 

L,l = [f,l] = jj^^n%, (47) 

where u, fi = x,y and summation over repeated indices is implied. This can be checked by 
inserting these equations into (Ja| • • Oo and comparing the result with the definition (p9D. 
(. . . I . . .)o is the Mori product with respect to Lq. Henceforth, bracketed indices refer to 
the magnetic field and unbracketed ones to the decomposition (^). Unfortunately, the 
conditions (0) are not satisfied in the Hubbard model. Instead, we derive with the help of 
Eqs. (|Tg)-(|g) and (0)-(|g) (see the appendix): 

[T, I] = ^-'"'^ < A, I] >o , (48) 



which should be equal to 



^Z^'y <l44l>o-J^X.. (49) 



This is obviously not the case. Similarly, we find 



C n ^ jT Sill Jx/rf. COS kr,, TIV' r ^ ^ -\ / \ 

T, J, = ^ f- < J,, Jy] >o , 50 

Ek. cos kx cos kyul^ 



instead of 



Efc^ sin kx n_ 



''ko 



Y^ka COS kx n' 



< [LJ^^ >o=if f^L- (51) 



ka 
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However, the conditions (|47| ) become exact in the continuum hmit or in the hmit of small 
band fillings. This is seen if we take explicitly into account the lattice spacing a in the 
arguments of the trigonometric functions which we tacitly have set equal to 1. Then we may 
expand cos k^a ~ 1 and sin/cj^a ~ k^a which proves the statement immediately. Thus, the 
violation of the conditions ( P7| ) on the lattice refiects its reduced symmetry in comparison 
with free space. Since the conditions (^) are properties of the free model, we may assume 
their approximate validity without taking the risk of not describing the local interactions 
(I) correctly. 



Before proceeding, we show how the conditions (|47|) appear within the formalism outlined 
in Ref. |T^. We expand Eq. (|3^) in the frequency regime, where the expression \Xxx{z)/x^\ 



is very small, i.e. for high enough frequencies and use a couple of times equations of motion 
for correlation functions A; B Thus we may show that M{z) can be represented as 
follows (cf. Eq. (H): 

2X°M(z) ~ ^ ^^^^ _ (52) 

Z 

The first term will be investigated in the next subsection and turns out to be regular for all 
frequencies; R{z) can be written as: 

R{z) = < 4; [T, Jy] - < [f , 4]; Jy >, +2nlyXxx{z) 

_2 <[X, Jy] >o 4>xx{z) - (pxxiO) 



(53) 



where 4>xx{z) =^ Kx] Calculating the function (f)xx{z) following the lines outlined in 

the next subsection, we may prove that R{z) is indeed divergent in the limit z 0, which, 
however, is an artefact of perturbation theory. The first representation of R{z) in the last 
set of equations shows that the condition (^) implies R{z) to vanish identically, if we take 
into account the symmetry properties fi^^ = —^y^ and Xxx{z) = Xyy{z). 
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B. Reduction to ordinary correlation functions 



We are interested in the memory function appearing in Eq. ([3^) , whose Laplace transform 
is given according to Eq. (|30D as 



M(t) = ^^{QLJ^\e''i'^^'\QLJy) . (54) 

Due to the approximation (|T7|), the free part Lq of the Liouville operator does not contribute 
to the operator QL\Ji,). Hence, to leading order in the interaction strength f/, we obtain 

M{t) = ^{QoLMe''^'''^''^''\QoLjy)o. (55) 
X 

Since Qo commutes with Lq and because of the idempotence of Qo [Jl^, we may free ourselves 



of all superprojectors Qo with the exception of one, say, that within the "ket" IQoLiJy). 
However, even this last appearance of Qq may be omitted, since its part Pq leads to a term 
proportional to the following first order expression of the frequency matrix = 
g^)^{Jx\Li\J^)o. Here summation over equal indices is implied. However, the frequency 
matrix is easily traced back to n^^ =< c^^ >, whose first order contribution vanishes. 
Thus, vanishes as stated and we arrive at M{z) = {P/x^){Kx\i/{z + Lo)\Ky)o, where 
we have defined 

K = [V, I] . (56) 

With the identity z/{z + Lo) = I—Lq/^z + Lq) and the symmetry property^ Kx] Ky ^°2= 
— ^ Kx', Ky which may be traced back to Eq. ( p7|) by means of two equations of motion, 
we eventually arrive at 

,/M(z) = - «-^-'-^^»° . (57) 



Now, we must evaluate this correlation function for the free tight binding model (0) to 
leading first order in the magnetic field. In order to derive explicit expressions for the 
operators (|56|), we introduce the following combination of Blochoperators: 
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. 4 4 ^ Cr, . (58) 

This is the basic building block for the operators Ky. To see this, we insert Eqs. (|^) and 
(10)- ([T^) into the definition ( ^61) and write the result in terms of Blochoperators. We find: 



bk, = Het^ ^ E [C{ki, h + q) - C{h - q, h)] (62) 

kik2kqcr 



K, = + 6Ky , (59) 
K^'^ = ^ [BAh. + q) - By{k, - q, k,)] (60) 

kik2 kqa 

Sk, = 0, (61) 

u 

where the matrices By{ki,k2) and C{ki,k2) are defined in the appendix, cf. Eqs. (|A2|) and 
(|X3|). Since SK^ vanishes, the expansion of the correlation function of Eq. (^) to first order 
in the magnetic field reads: 

« k,, ky »° = C'\z) + C'"{z) , (63) 
C^^(^) =<E:f >f) , (64) 
C^^^(^)=«E:f;Kf»f) . (65) 

Obviously, it is sufficient to calculate the correlation function consisting of operators ( |SHD 
within the tight binding model (|T^, however, to first order in the magnetic field. But first, 
we note that the functions (^) and ( p^D are two- and three-point correlation functions, 
respectively. This is explicitly seen within the Matsubara representation where the expansion 



of ^ kx] ky up to first order in the "perturbation" (|l^) yields: 

1 ri^ 



C'\iu^) = -- drdr' < Tr {i^f (r) 6ky{r')} >[,°) e^-'"(--') , (66) 



C''\iu^) = ^j^j^j^drdr'dT" <T,[k^^\r)kf\T')5f{r")}>f e^-'"(--') . (57) 



C. Expansion to first order in the magnetic field 



As already mentioned, our next goal is to calculate the correlation function generated 
by the operators ( ^Sf ) up to first order in the magnetic field. This is accomplished by means 
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of its equation of motion with respect to the tight binding hamiltonian (|T3]): 



Here and in the following, the index is omitted. The expansion of the expectation value 



on the Rhs. with respect to the magnetic field is standard [|17l and yields 



^ ^ ["^fcl.telfclq ' k[,k'2\k'\q'\ >^ ~ < ["^fci,fc2|fc|g ' fci.fc^jlfe'lg*'] ^ ' '^^^^ 



where v arises from the expansion of the ^'-matrix related to the "perturbation" (|T5D and is 
therefore given by 

v= dr e-(^'"'-'^^) 5f e-<*'''-^^^ . (70) 
Jo 

With the representation of 5T in terms of Blochoperators, 

5f = -zetH X: D{h, hYl^cj:^, , (71) 

k\k2C! 

where the matrix D{ki,k2) is also defined and further evaluated in the appendix (cf. Eq. 
(|A^)), we find more explicitly: 

^ ^ _ 1 

V = -tetH Y: Dih, h) _ cl^c^,^^ (72) 

Inserting the expansion (^9]) into Eq. (|68D , we obtain the following zeroth and first order 
terms for the correlation function to be determined: 



« ^ ; 4' P ,P, = '^''^"''^ , (73) 

ki,k2\k\q ^ k[,k!2\k'\q' ^ r -|- f ^ ^ — f ^ -|- f ^ — ^ ' 



k-q ^k I '-fci '-fc2 



fei,fc2|fe|g ' /c'fcilfe'k' ^ r -I- _ -I- _ f ^ 

^ ''k-q ^k ^ ^ki *=fc2 



(74) 
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The second term on the Rhs. of Eq. (^) still contains a correlation function. Fortunately, 
this function is related to the hamiltonian (0) without magnetic field thus being directly 
reducible to expectation values by means of its equation of motion: 

« l^T, ; »f ^ ^ J'^'Z :'f^ • (75) 

^ ^ ^ ^k'-q' + ^k' ^k\ + ^k' 



In summary, the problem of calculating the memory function ( ^Tf) to leading order in the 
magnetic field has been reduced to the calculation of expectation values within a free tight 
binding model without magnetic field: The relevant information is contained in the equations 



(pO|), (§3)-(^) and ([731)- ([75|) . The rather cumbersome calculations are roughly sketched 
out in the appendix. We write the memory function contribution to the Hall constant (cf. 
Eq. (ID) as follows: 

e \ 2ta2 



where a„ is the amplitude of a nearest neighbour hop and is related to the static susceptibility 
(PB|) via x° = 4:tNan. This in turn implies: 

= ^Ecos/co^rig,. (77) 

fecr 

m^\z) and m^^\z) arise from the two- and three-point correlation functions, respectively 
and are represented with regard to the further strategy as energy integrals: 



m'\z) = -^J *i / de2 J de[J de', J(^|ei, esje'^, e^)L^^(ei, esle'i, e'2) , (78) 

Vd^ 

Q{z) = J de J dei j de2 J de[ J de 



m'^\z) = ^(Q{z)+Q{-z)), (79) 



, ^(2:|e, e2|ei, 4) - ^(^l^i, e2\e[, e'^ 



ei 



X L^^^(e|ei,e2|e;,e^) . (80) 

The integrands feature the following abbreviations: 

. ^ _ /(ei)/(e2) [1 - /(e;)] [1 - m)] - [1 - /(eQ] [1 - /(e^)] /(e;)/(e- ) 



{cos + COS - [sin k^, - sin k[^] V cot( ^^" ^ ~ ~ ^^" ) } 
{cos kly + cos k2y — cos k'^y — cos /cg^} 

27r5(A;i, + A;^, - k[y - k'^^) (82] 
X {cos kix + cos k2x + cos A;']^^, + cos /cg^; 

- [sin fc,. + sin k,^ - sin k[^ - sin A:^] P cot(^l^±^^^^^^-^)} 
X {sin kly + sin k2y — sin /e^^^ — sin k2y} 
X sinA;ij,27r5(A;iy + k2y - k[y - k'^y) ^ 



X 
X 
X 



;83) 



and /(e) = l/(exp(e — /i) + 1) is the Fermi function. < . . . >^ denotes the average over the 
first Brillouin zone, i.e. / df'k/ {2ttY{. . .). Note, that the last two equations ( ^21) and (|83D 
reflect the gauge fixation (P): In the ^/-direction, crystal momentum is conserved which is 
ensured by the 5-functions while the x-components of the Bloch vectors are coupled more 
complicatedly. In principle, we could do the momentum integrations numerically for a two- 
dimensional lattice and for given sets of external parameters temperature T, doping 6 and 
frequency uj. However, we may carry on our analysis a little bit by invoking a limit pioneered 
by Metzner and VoUhardt in the context of strongly correlated electrons [0, namely the 
limit of infinite spatial dimensions. In this limit, the momentum integrals decouple and we 
are left with energy integrals over smooth functions. This procedure will be discussed in the 
next subsection. 



D. The limit of infinite lattice dimensions 

We may question the relevance of this limit, since the important physics of the high- 
Tc superconductors is known to take place in Cu-0 planes. Many authors have addressed 
this issue and much evidence has been revealed in favour of the relevance of this limiting 
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procedure even for two-dimensional systems, see e.g. Ref. |T9|. Instead of immersing ourselves 
in this debate, we take the following point of view: The main reasons behind the anomalous 
properties of the high-Tc materials seem to be, firstly, the strong electronic correlations and, 
secondly, the two-dimensionality of the relevant Cu-0 planes. Taking the limit d —>■ oo helps 
us to separate the impact of the correlations and to suppress effects of low- dimensionality 
like e.g. van Hove singularities. In this sense, the limit ^ oo is interesting in itself. 

For our problem, the most important aspect of the limit (i — ^> oo is the following: For the 
Hubbard model to retain its nontrivial dynamics, the parameter t has to be scaled properly 
with d according to 

(in this subsection, we set t* = 1). Only then does the Hubbard model capture simul- 
taneously the itinerant and local aspect introduced by the hopping and interaction term. 



respectively. On the other hand, we are tempted to conclude from the scaling (|8^ that any 
transport stops to be possible in c? = oo. In fact, a more thorough investigation shows that 
the longitudinal and the Hall conductivity are of order 1/d and l/rf^, respectively. But this, 
in turn, implies that the Hall constant remains finite in d ^ oo. In the following, all we 
need to know is how to calculate averages over Brillouin zones of the type ( p2D and (pSj). 
The corresponding procedure is explained in the appendix and enables us also to calculate 
simpler quantities as the density of states D{e) of the band model (0), the nearest neigbor 



hopping amplitude ( ffTl) and the amplitude of a hop diagonally across the unit cell, i.e. 

«rf = ^Ecosfcxcos/cj^ng^. (85) 

ka 

In the case of the amplitudes (^) and (|85D , the functions 

A(e) = < cosA;,.5(e - e,^) >,^= -^^(e) , (86) 
5(e) = < cos K cos ky6{e - e^) >g= ^(e^ - i)D(e) (87) 

come into play. Essentially, they are Gaussians multiplied with the first and second Hermi- 
tian polynomial, respectively, since they are derivatives of the Gaussian density of states: 
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D{e) = ^e-^\ (88) 



The functions (^2]) and (|83D are found to be combinations of the functions (|8q) and (^3) 
may be written as 

L^^(ei, e2|e;, e^) = Z}(ei)Z}(e2)D(e;)Z}(e'2){e? + eiea + ese'i - - - e^ei} , (89) 
L^^^(e|ei, eale'i, e^) = D{e)D{e^)D{e2)D{e[)D{e'^){-e + ei + + e'^ + e^} . (90) 

In order to handle the singularity of the Hall constant in the empty band limit correctly, 
we shall discuss the perturbation-theoretical results for the Hall constant normalized to its 
f/ ^ limit: 

?^1^2E1 = 1 + u\K^ + K'\z) + K''\z)) . (91) 
Here, the t/ = Hall constant is given by 

= i (^ 

(the subscripts indicate [/ = as above) and K^\z) and K^^^{z) arise from the func- 
tions ( |78D and ([79|), e.g. K^\z) = {2d/ ado)m^\z). K°° is the perturbation-theoretical 
contribution of the infinite frequency Hall constant (pS]). Since the latter is given by 
RJI = (l/e)(2arf/(2a„)^), it suffices to calculate the density to second order in U. The 
derivation is standard and will therefore not be given here. Due to symmetries of the ex- 
pression (|90D , the function (^) may be simplified by means of various redefinitions of the 
integration variables: 



Q{z) = JdeiJ dt2j dt\j D{t{)D{t2)D{t[)D{t'^) 



X /(z|ei, e^le;, e',) {e^ + e', + e'J / de ^^''^"^ . (93) 



The terms of Eq. ( |9TD may now be evaluated numerically in the limit u —>■ on the basis of 
Eqs. (|7|), (|7|) and (0) along with the definition (^ and the results and (^). But first, 
we check whether the Hall constant (|9l| ) reduces to the familiar expression 1/ne, provided 
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the electron density n = {^/N)^^^n^^ is very low. We concentrate on zero temperature, 
where we find the Fermi energy to be given by = 1 + erf(ei7'). This implies ep —>■ —oo in 
the empty band limit. Then, the corrections on the Rhs. of Eq. (pl| ) vanish and the Hall 
constant is given by Eq. (p^). With 2a„o = D{ep)/\/d and 2ado = —eFD{eF)/d, we find in 
fact for n — 0: 

Rh = —^,^-. (94) 

eD{eF) en 



E. Numerical results 

First of all, we discuss the relative importance of the terms appearing on the Rhs. of Eq. 
(pip. Fig. H shows their doping dependence at T = (dashed lines) and that of their sum 
(solid line), where the doping parameter is defined as 5 = 1 —n. All functions vanish in the 
empty band limit (6-^1) and exhibit monotonic behaviour with decreasing doping. This 
refiects the fact, that the suppression of doubly occupied sites introduced by the Hubbard 
repulsion (^) becomes more effective with increasing electron density. As for the signs of 
the three contributions, only that of the two-point correlation function is negative. This, 
however, is sufficient to render the sum of all terms negative (solid line in Fig. |l|). This 
remains valid at finite temperatures: The term K°° is positive for all temperatures and 
dopings considered, but is always overcompensated by the memory function contribution 
K'\iO^) + K^^\iO~^). Thus, our perturbation-theoretical results clearly indicate the ten- 
dency of the Hall constant to change its sign below half filling. However, for this to happen, 
the memory function contribution must be taken into account. 

To study the doping and temperature dependence of this precursor effect in greater 
detail, we shall extrapolate Eq. (pl| ) to correlation strengths U big enough for the Hall 
constant to exhibit a sign change. Ultimately, we fix U such that this sign change occurs in 
the parameter regime observed experimentally in the case of the compound La2-5Sr5Cu04. 
We shall measure U in terms of the bare band width W, which may be chosen, due to Eq. 
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(^), as W = 2t* (from now on, the hopping parameter t* is exphcitly taken into account). 
Then, it follows from Eq. (^) and Fig. |1|, that the T = threshold U, above which the Hall 
constant becomes positive, is U = 2.181^ at half filling {6 = 0) and increases monotonically 
with increasing doping and ultimately diverges in the empty band limit 6^1. 

Before proceeding, we touch upon the issue of how to relate our theoretical results to 
experimental measurements. Firstly, the d = oo hopping parameter t* has been estimated 
crudely as 0.2t* ~ 500K ||2^. Secondly, we shall express the Hall constant in units that allow 
direct comparison with experimental results. This requires that charge carrier densities are 
taken with respect to the volume of a unit cell. On the other hand, the electron density 
n, appearing in Eq. (0), denotes the average number of electrons per lattice site. From a 
theoretical point of view, this definition is convenient since it is independent of the lattice 
dimension d and remains meaningful in the limit d oo. Therefore, in order to compare 
our theoretical results with measurements on a certain cuprate, we have to multiply the Hall 
constant of Eq. ( pi]) by Q/u, where fl is the volume of a unit cell and u the number of Cu 
ions therein. In the case of La2-5Sr5Cu04, Q = 186 and z/ = 2. 

After these preliminary remarks, we proceed with the investigation of the intermediate 
correlation regime, where a sign change is possible. Fig. ^ shows the doping dependent 
Hall constant for several temperatures and the choice U = 2.3W (solid lines) as well as two 
experimental curves for polycrystalline samples of La2_5Sr5Cu04 taken from Ref. ^ (in- 
set). We see that the sign change occurs close to 5 ~ 0.3 for temperatures below 300K, in 
agreement with experiment. For lower dopings, our Hall constant exhibits a maximum and 
ultimately vanishes at half filling, irrespective of temperature. This reflects the fact that, 
in our perturbation-theoretical result (|9T|), the Hall constant of the bare band is merely 
renormalized by a finite factor. Such a factor may change an overall sign but never can turn 
a vanishing quantity into a nonzero one. Thus, our results fail to account for the observed 
1/6 law of the Hall constant near half filling. And that is why our perturbation-theoretical 
curves cross the zero line (dotted line in Fig. ^ with slopes that are two orders of magnitude 
smaller than those of the experimental curves. Apart from this deficiency of perturbation 
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theory, the other dependences are in quahtative agreement with experiment. In Fig. ^ the 
temperature dependence of the Hall constant is shown for various dopings within the range 
0.1 < S < 0.4 (solid lines) and compared to experimental results from Ref. [|11|], again for 
polycrystalline samples of La2-iSr5Cu04 (dashed lines in the inset). Despite the already 
mentioned difference in the order of magnitude of the Hall constant, our curves display the 
same features as the experimental ones: A maximum at low temperatures followed by a 
regime in which the Hall constant decreases monotonically up to unusually high temper- 
atures. We have not been able to determine the exact location of the maximum due to 
numerical difficulties at nonzero temperatures below 0.05t*. Experimentally, it occurs above 
Tc. The fact that it appears within the Hubbard model suggests that it is not related to the 
onset of superconducting correlations. This is further supported by comparing the data of 
"90-K" and "60-K" YBa2Cu306+x [§]■ ^^is compound, the location of the maximum in 
Rh(T) does not depend on the doping values corresponding to the range 60K < Tc < 90K. 
As for the decrease of the Hall constant as a function of temperature, it is experimentally 
found to be most pronounced at optimal doping, i.e. 5 = 0.15 in the case of La2_5Sr5Cu04. 
Our results exaggerate the doping range where this decrease is markedly visible. At least. 
Fig. ^ shows that the decrease is least pronounced for the curve corresponding to the lowest 
doping value, 6 = 0.1. Furthermore, at high dopings, where the memory function contribu- 
tion becomes unimportant, the Hall constant becomes almost temperature independent. 

What about the observed quadratic dependence of the Hall angle on temperature for 
small dopings |^,0,|1^? This law cannot be verified on the basis of Eq. (|9T|) alone, although, 
it cannot be falsified either. To make a check on this law, the longitudinal conductivity is 
needed as well. In principle, the calculation of this quantity can be done along the same 
lines leading to Eq. ( |9lD and is left for future work. 
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IV. CONCLUSIONS 



In summary, we have devised a memory function approach to the Hall constant in strongly 
correlated electron systems, which enables us to cope directly with the Hall resistivity. As 
a first step, its usefulness was demonstrated in a perturbation-theoretical treatment within 
the single-band Hubbard model. To obtain a regular expression for the memory function 
contribution for all frequencies, we assumed that the subspace of the operator space spanned 
by the current operators Jx and Jy is invariant with respect to actions of the unperturbed 
Liouville operator. This approximation was shown to become exact in the continuum limit. 
Furthermore, it affects only the properties of the unperturbed system, i.e. the tight-binding 
electrons. Therefore, we do not expect the omitted terms to change the physics in an essential 
way. On the basis of this approximation, the memory function was calculated to leading 
order in the correlation strength U and shown to decompose into a two- and three-point 
correlation function, when expanded to first order in the magnetic field. The complicated 
expressions obtained for these two functions were simplified considerably by invoking the 
limit of infinite spatial dimensions. While this approximation still catches the impact of the 
correlations, it smoothes out effects of low dimensionality. Except for the doping dependence 
of the Hall constant in the vicinity of half filling, we have been able to reproduce the unusual 
experimental findings in connection with high-Tc superconductors as La2_5Sr5Cu04. In 
particular, we could explain the sign change of the Hall constant as a function of doping 
and its decrease as a function of temperature up to unusually high temperatures. However, 
we had to chose U = 2.3W {W: bare band width), which is, strictly speaking, outside the 
perturbative regime. Since the cuprates are believed to undergo a transition from a Fermi 
liquid to a strong-correlation regime when the doping approaches its optimal value from 
the overdoped side 0, it is not astonishing that the 1/6 law in the vicinity of the mother 
compound cannot be described within perturbation theory. A treatment of the Hall effect 
in the opposite limit of strong correlations is in progress [jl2[ . 
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APPENDIX A: LATTICE SUMS OVER NEAREST NEIGHBOURS 

Since our gauge choice, cf. Eq. (|^), permits periodic boundary conditions in the y- 
direction only, we must carry out all sums over nearest neighbours according to the following 
formula: 

^ HiR,R + 5)=Y.Y.T.H{R,R + Syy)+ E T.HiR,R + S.£) ■ (Al) 

R,R+S Rx+6x Ry 

In the first term on the Rhs., we may carry out the sum over Ry and by independently. As 
for the second term on the Rhs., we must make sure that R^ and R^ + are neighbouring 
elements of the set {i?™", + 1, • • • , RT""} with i?^^"^' = i?™" + A^"^ - 1 as explained in 
the text following Eq. (p!2D. 

First of all, the matrices appearing in Eqs. (|60|), (|6^ ) and ( [7TD are defined as follows: 



B,{k^M) = % E Ke-^^^^^^'^^^^-^ , (A2) 

C{kM = \ E (5,)'i?.e-^^^^(^+^>^^^^^^ , (A3) 

D{k^M) = ^ 5,i?,e-^^^^(^+^>^^^^^ . (A4) 

In all these cases, the integrand is proportional to a component of a nearest neighbour 
vector, which entails that one of the two terms of Eq. ( |A1| ) vanishes. The summation over 
the ^/-components Ry and by is straightforward. In the case of Eqs. (|A3| ) and (|A4[), the 
following sum is to be calculated: 

1 pi{k2x-k^x)RT'^ 

-Y.R. e^^'^--'^-)^^ = b,^^\,,Rl + (1 - .(,,^_,^^)_^ , (A5) 
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where R^. = + -R™"^)/2 is the x-component of the lattice's center of gravity. In 

summary, we obtain the following results: 

I — 



D{ki,k2) = -2ismkiy6k,^ik2^ ih.^lk^.K + (1 - '^fci.lfcaJ ^^(fc.^.fc,^) _ J ■ (A9) 

In the following section, we shall see that, once these expressions are inserted into observable 
quantities, the final results become independent of the lattice location. 
The commutators ( ^8|) and (|50|) are also calculated using Eq. (|A1|) . 

APPENDIX B: EVALUATION OF THE TWO- AND THREE-POINT 

CORRELATION FUNCTION 



We begin by inserting Eqs. (BUI) and (|63) into Eqs. (01) and (pSf) and by using Eqs. (|73| 



(fTSP along with Eqs. (|A^)-( |A9| ). Then, the correlation functions ( [6^ ) and (|65| ) are given as 
sums over terms that contain expectation values with respect to the momentum conserving 



Hamiltonian ([1^) . Taking into account the corresponding delta functions, we see that the 
diagonal elements of the matrices (|A6|) , (|A8|) and ([A9|) lead to vanishing contributions due 
to symmetry arguments. Moreover, the combination of all exponentials whose arguments 
are proportional to may be replaced by one due to the delta functions, expressing 

momentum conservation in the model (p!4D. With the help of the function 

1 



S{k) = Ny6k^\o{l-h.\o) 



I . r „ ,k. 



^Ny6,^\,[l + tV coi{^)] (Bl) 



and after some straightforward manipulations, we arrive at: 

-2iHet'^NU'^ = jyl hiylk^y-qy S{k'2 + q' - k{] 

k\k2kqa k'^k'^k'q'a-' 
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X 



L fci,fc2|fe|ij ' k[,k!2\k'\q' 



.(0) 
^0 



^lW = ]^ E E E •5(P^2-Pl)54+-,.fc||g 

A^lfc2fc(^o- k[k'^k'q'u' PiP2T 

■ , [sin fc;^ - sm{k[y - q'y)] (1 - e"*^-) (e'^^i- + e-^^^^) 



X 



sin Ply 



X 



g/3(epi-<:p2) — 1 



'1 ^P2 ' ^fe- 

A(^) = ]^ E E E^bi 

klk'ikqa k'^k'^k'q'a' P^P'^'^ 



'PIT 

~ Pi) ^k'^+q'~k[\(i 



^ sin Ply [sin k[y - sm{k[y - q'y 

{^+^k-q-^k+^k,-%) 

^0 



^ ^k'-q"' 



jihix _|_ g" 

^k' ~ ''k[ 



-ik2x 



g;jJ(l-e-^'?-)(e^^ 



Pl|fc2 



^0 



+ Sr\-o 



^P2\k—( 



\- kxM\k\k-Vx ' k'^,k'^\k'\q'\ ^0 



"pilfc ^ I- fci,fc2|P2|P2-fc+-f ' "^^.yfc^lfc'lg'-l ^0 



From the definition f^), we find: 



(B2) 
(B3) 



(B4) 



(B5) 



The further evaluation of Eqs. (|E^ ), (B4) and (B5) requires the calculation of expecta- 
tion values. Fortunately, not all terms that arise from the corresponding factorizations, 
contribute: Terms proportional to or ^i^,^^ may be omitted in any case. And terms pro- 

28 



portional to 5j^,^^,_^, and do not contribute in the case of Eq. (p^ ) and Eq. ( p4D along 
with Eq. ( |B5[ ) , respectively, due to the fact that S'(O) = 0. Thus, we see for example that 
a factor /(ep*J(l — /(cp-j)) may be split off from the expectation value within the integrand 
of Eq. (P^). This factor combines with the quotient that stems from the operator ( |72D 
according to 

The further calculation is not simple and takes some time, especially in the case of the 
functions (p^ ) and (|B5| ) . In the thermodynamic limit, where we may replace, e.g., Ny5ky\o ~^ 
2TT6{ky), {l/N)J2k ^ Ii Bz d-'^^/i'^'^Y etc., our final result may be written in the form of 
Eqs. m-<M)- 



APPENDIX C: BRILLOUIN ZONE AVERAGES IN INFINITE DIMENSIONS 



The calculation of Brillouin zone averages to leading order in l/d follows the procedure 



outlined in Ref. |21]. In the following, Schlafii's integral representation of the Bessel functions 



will play an important role: 

< e^-^+*^™^'= >k= J|„|(r) = Gn{r) . (CI) 

Here, n is an integer and we have used the notation < . . . >k= /^^ ■ ■ • dk/{2'K). By means 
of the Fourier representation of the delta function, we may write, e.g., the function defined 
in Eq. (|8^) as follows: 

/OO fj c 
_^l^V-Hr)f[Ur)r'e-\ (C2) 

where here and in the following, r = s/\fd. Expanding the Bessel functions in powers of r 
and taking the limit d ^ oo, we find: 

B(.) ^ . (C3) 
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Thereby, we used -D(e) = J^ooT' ^ "4"*"*^^, which is derived analogously. This proves Eq. 



(0). The corresponding evaluation of Eqs. ( ^21) and (^) requires a Fourier series expansion 
of the cotangent: 

Vcot{^) = -Y.sgniR)e'''' . (C4) 

Here, the sum is over all integers R, except for the zero. Since we are working in the 
thermodynamic limit, R may be considered to be a component of a lattice vector. To prove 
this representation, we start out with a known formula for the coefficients 6„ appearing in 
the ansatz Pcot(|) = Y.T=i sin(n/c): 

| = ^£§cot(|)sin(nfc). (C5) 

With the substitution z = e'*^, we may perform the principal value integration by invoking 
the theorem of residues: The integration contour goes around the unit circle with the point 
z = 1 being excluded. Thus, we have to add the residue at z = to the half residue at 
z = 1. We obtain 6„ = 2 for all positive integers n, which proves the statement ( |(J4D . In the 
following, we show how the Rhs. of Eq. (^) is evaluated to leading order in 1/^/d. Writing 
the delta functions, that contain energies, in terms of Fourier integrals and introducing an 
additional momentum average < (27r)'^5(fc + ^2 — k'l — k'2) . . . >^ = 1, the problem reduces 
to the calculation of the following expression: 



d 

x2n6{k^ + k2x - - k^^) 2TT5{ky - kiy) W 27i6{kj + k2j - k[ - - k'^j) 

i=3 

X {cos kix + cos k2x + cos k[^ + cos k2^ — [sin ki^ + sin k2x — sin k[^ — sin ^2^] V cot( — ^—^ — -)} 
X {sin kiy + sin k2y - sin k[y - sin fc^ J sin kiy2n 6 {kiy + k2y - k[y - /c^^) )kkAk[% . (^6) 

Since every dimension j contributes a term — cosfcj to the band dispersion e{k), 

the expression ( |C6| ) decomposes into d factors. For example, one factor arises from all 
x-components: 
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^gi[rcosA:.+ricosfci.+r2Cosfc2.+ricosA:;^+r^cosfcy 27l5{k^ + - - k'^^) 

X {cos kix + cos k2x + COS k[^ + cos k2x 

-[sin /ci^ + sin fcsx - sin k[^ - sin /c^^] P cot(— ^^-^)} )k^k^,k2.k[^k'^^ . (C7) 

This expression decomposes further into eight terms corresponding to the ones in the curled 
brackets. Each has to be evaluated with the Fourier series expansion ( P^ and that of the 
delta function: 

oo 
R=-oo 

For example, the term — sin ki^ V cot ( ) gives rise to the contribution 

^ ^sgnHG'„_„(r)i[G'„+i(ri) - G'„_i(ri)]G'„(r2)G'^K)G'„(r^) , (C9) 

where Eq. ( |CT1) has been used. The leading order in 1/ \/d reads: 

- Gi(r)Go(ri)Go(r2)GoK)Go(r^) . (CIO) 

This contribution to the expression (|C^ ) is of order l/y/d, as are all the others. The sub- 
sequent Fourier integrals over the variables Sj yield combinations of the functions (|86|), (|87D 
and (|88|) thus leading ultimately to Eq. (BOf). Eq. (R^) is proven analogously. 
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FIG. 1. The corrections of Eq. (M) at T = 0. The curves A, B, C, and D represent the functions 
K^^{5), K^^\5), and their sum, respectively. 




FIG. 2. Hall constant as a function of doping for U = 2.3 IF. Inset: Data for polycrystalline 
samples of La2-5Sr5Cu04 taken from Ref. ||5|]. 



34 




35 



